Efficient high-fidelity quantum computation using matter qubits and hnear optics 
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We propose a practical, scalable, and efficient scheme for quantum computation using spatially 
separated matter qubits and single photon interference effects. The qubit systems can be NV-centers 
in diamond, Pauli-blockade quantum dots with an excess electron or trapped ions with optical 
transitions, which are each placed in a cavity and subsequently entangled using a double-heralded 
single-photon detection scheme. The fidelity of the resulting entanglement is extremely robust 
against the most important errors such as detector loss, spontaneous emission, and mismatch of 
cavity parameters. We demonstrate how this entangling operation can be used to efficiently generate 
cluster states of many qubits, which, together with single qubit operations and readout, can be used 
to implement universal quantum computation. Existing experimental parameters indicate that high 
fidelity clusters can be generated with a moderate constant overhead. 

PACS numbers: 32.80.-t, 78.70.-g 



Quantum computation (QC) offers a potentially expo- 
nential computational speed-up over classical computers, 
and many physical implementations have been proposed. 
Particularly promising proposals are those in which uni- 
tary operations and readout in matter qubits are imple- 
mented via laser-driven optical transitions. Examples 
are the original ion-trap proposal JjJ , NV-centers in dia- 
mond 2] , and schemes utilizing the Pauli-blockade effect 
in quantum dots with a single excess electron 'i', . Sin- 
gle qubit operations and readout, using a combination of 
optical and RF control fields, have already been demon- 
strated in ion trap and NV-diamond systems IM l^i 
while a number of promising techniques for optically ad- 
dressing quantum dot spin qubits have been proposed 
I^Q. In all these cases, the ratio of the single qubit op- 
eration time to the intrinsic decoherence times suggests 
that very high fidelity operations are possible. 

However, there are substantial difficulties in scaling 
these implementations to the large numbers of qubits re- 
quired for useful QC. Multi-qubit gates are facilitated by 
a direct interaction between qubits. Thus adding a new 
qubit to a quantum register, together with the associ- 
ated control fields, necessarily modifies the Hamiltonian 
of the system. This can mean that, as more qubits are 
added, logic gate implementations become progressively 
more complex, and furthermore, new decoherence chan- 
nels can be introduced. Furthermore, the need to opti- 
cally address individual qubits (e.g. in NV diamond or 
quantum dot systems) can lead to seemingly contradic- 
tory system requirements: the qubits need to be suffi- 
ciently well separated to be resolved by the optical field, 
but must be close enough such that two-qubit logic can 
be implemented via the inter-qubit interaction. 

A potentially promising solution to these scaling chal- 
lenges is to perform distributed qunatum computing, in 
which the matter qubits are spatially separated. In this 
case, there is no direct interaction between the qubits. In- 
stead, Entangling Operations (EOs) between qubits are 
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FIG. I: In the circle: The qubit system {|T),|i>} with the 
excited state |e). The 7r-pulse affects only the transition |J,) 
|e), and the emission of a photon into the cavity mode brings 
the excited state back to the qubit state Using a pair of 
leaky cavities with this system and conditioning on a single- 
photon in the detectors D+ or D- creates entangled qubit 
states. The 50:50 beam splitter (BS) is used to erase the 
which-path information. 



implemented via single photon interference effects. A 
number of schemes to entanglepairs of distant qubits in 
this way have been proposed [1 S H H [H [l2 . Re- 
cently, it has been shown that unitary logic gates can 
be can be performed in this manner |l3L Il4|. However, 
the latter schemes are either inherently non-deterministic 
jl^ or are sensitive to photon loss or photodetector in- 
efficiency |l3j and it is not clear whether they can be 
used for scalable QC. Other schemes using single photon 
interference effects together with local two-qubit unitary 
operations have also been proposed [21 Eil ■ 

In this Letter, we propose a fully scalable scheme for 
distributed QC using individual matter qubits assuming 
only single qubit operations. Our scheme is robust to 
photon loss and other sources of errors, and uses opti- 
cal transitions of the qubit system, together with linear 
optics and photodetection to entangle pairs of spatially 
seperated matter qubits in a non-deterministic manner. 
A key observation is that even such a non-deterministic 
EO is sufficient for scalable QC: our EO can be used to 
efficiently generate cluster states of many qubits, which. 
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together with single qubit operations and measurements, 
are capable of universal QC the context of linear 

optics QC 0], it has recently been shown 0, 0| 
that the cluster state model can be used to significantly 
reduce the resource overheads required for scalable QC. 

We consider matter systems comprised of two long- 
lived, low lying states |t) and ||), and one excited state 
|e), in an L-configuration (see Fig.P). The system is con- 
structed in such a way that an optical 7r-pulse will induce 
the transformation ||) |e) and |t) |T)- The transi- 
tion It) ^ |e) is forbidden, e.g., by a selection rule. The 
states It) and ||) represent the logical qubit states |0) 
and |1) respectively. We assume that high fidelity single 
qubit operations and measurements can be performed on 
these logical qubits. Physical systems that have a suit- 
able level structure include NV-centers in diamond 0, 
quantum dots with a single excess electron 0, 0, and 
various trapped ion and atomic systems. Each such sys- 
tem is embedded in a separate optical cavity, such that 
only the ||) ^ |e) transition is coupled to the cavity 
mode. One end of each cavity is leaky, with the leakage 
rate of the i**^ cavity given by 2^;. The light escaping 
from the cavities is mixed on a 50:50 beam splitter, the 
output modes of which are monitored by two vacuum- 
discriminating detectors, -D+ and with efficiency rj. 

The scheme proceeds as follows. Firstly, both qubits 
are prepared in the state |-|-) = (|t) + |i)) /\/2 using lo- 
cal unitaries. We then implement the following sequence 
of operations: i) Apply an optical 7r-pulse to each qubit, 
coherently pumping the population in the ||) state into 
the |e) state; ii) Wait for up to a time twait for a photo- 
detection event in either or in) Wait for a fur- 
ther time trciax for any remaining excitation in the qubit- 
cavity systems to relax; iv) Apply an X operation to both 
qubits, coherently flipping the spins as |t) \l) and 
II) It); v) Repeat steps i-iii. 

Appropriate values for twait and treiax are determined 
by the system parameters discussed below. If zero or 
two photo-detection events are observed on either round 
of the procedure, the scheme failed, and the qubits must 
be newly prepared before re-attempting the entangling 
procedure. On the other hand, if one (and only one) 
photo-detection event is observed on each round of the 
protocol, the scheme has succeeded, and a maximally 
entangled state is prepared with unit fidelity (given ideal 
systems). We call this technique double heralding, and it 
turns out to be exceedingly robust against most common 
experimental errors. 

We analyzed the scheme in detail using the quantum 
trajectories formalism For clarity, we first consider 
the ideal case, in which the detectors have unit efficiency 
(77 = 1), and spontaneous emission of photons from the 
transition |e) ||) into modes other than the cavity 
mode is neglected. During time periods where no de- 
tector clicks are observed, the conditional state of the 
system, in the interaction picture, evolves smoothly ac- 



cording to the effective Hamiltonian {h= 1) 

E |(U)..(e|cl+H. c.)-* J2 ^4S.- (1) 
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Here, gi denotes the Jaynes-Cummings coupling between 
the \e)i ^ transition and the mode of the i^^ cavity, 
and Ci is the corresponding annihilation operator. For 
the purpose of illustrating the ideal case, we assume that 
systems A and B are identical, such that gA = 9b =9 
and ka ^ Kb — K, and that k > g. 

When a single click is observed in detector D±, the 
state of the whole system discontinuously evolves as 
IV'(t)) c±\tp{t)), where c± = {ca ± Cb)/V^ denotes 
the corresponding jump operators. Thus, after steps i- 
ii of the entangling protocol, conditioned on observing a 
detector click at time ti < iwait, the unnormalized state 
of the whole system is 



\i>{ti)) = a{ti)\^±)+a{ti)l3{ti) 
+2a^{ti) 



U,0;e,0)±|e,0;i,0) 



2, Ji,0;i,l)±U,l;t,0) 



^/2 



V2 



(2) 



Here, \qA,PA', IBiPb) is the state of the whole system, 
with qA(B) and Pa(b) denoting the states of matter sys- 
tem A [B) and cavity mode A [B) respectively, \^^) = 
(||,0;t,0) ± It, 0; t, 0))/\/2 are maximally entangled 
states, a{t) = -ig/{AyjK? - g'^) (e-r=iowt/2 _ e-rfastt/2) 
and /3(^) = ^(1 + K/yZ/t^ -g2)e-r3iowt /2 ^ i( i _ 
k/ _ (^2^g-rfaEitt/2^ where Ffast = k -f \/ — g"^ and 
Tsiow — \/ ^ ■ In order to obtain a detector click 
with significant probability, iwait should be chosen to be 
a few times V-^l^. 

Equation (|2Jl implies that it may be possible to observe 
a second detector click on the first round of the pro- 
tocol. However, realistic photo-detectors typically can- 
not resolve two photons arriving in quick succession [23| . 
Within the quantum trajectories description, this can be 
simulated by assuming that no information is available 
from either detector after the first click. After a time 
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^ r , , the system decoheres to the state 
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where = l + |/3(ti)p-|-|a(ti)p. The undesirable second 
term in Eq. is removed by applying steps iv and v of 
the entangling procedure. If a photo-detection occurs on 
the second round, the final state of the system is a pure, 
maximally entangled state. If the two clicks are observed 
in the same (different) detector(s), the final state is |^+) 
(1 4'")). Each of the four possible successful outcomes 
occurs with probability 1/8, leading to a total success 
probability oi p = 5. 

We also analyzed the scheme in the non-ideal case, al- 
lowing for less than perfect detector efficiency ry < 1, and 
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finite spontaneous emission into free space (71 = 72 = 
7 > 0). These imperfections do not reduce the fideUty of 
the final state, but do reduce the success probabiHty (see 
Fig. Note p has a quadratic dependence on 77, while 
p decreases rapidly for 7 > Fsiow 

The dominating experimental imperfections that do 
reduce the fidelity can be classified into three groups: 
(1) decoherence of the matter qubits; (2) dark counts in 
the detectors; and (3) imperfect mode matching of the 
photons incident on the beam splitter. Firstly, the ef- 
fect of spin decoherence depends on the way the cluster 
states are generated, and can be estimated by compar- 
ing the spin decoherence time with the "clock time" 
tc ^ lOr^jiJ^ at which the EO can be repeated. If the 
preparation of cluster states is performed in parallel, the 
typical time overhead is m clock cycles (see below). Thus 
the average age of a qubit the moment it is added to the 
cluster is (m/2)ic, and to < 8 for reasonable detector effi- 
ciencies. Assuming a reasonably strong cavity qubit cou- 
pling, g — IOO7, and critically damped cavities {g ~ k), 
the size of errors due to spin decoherence is given by 
e ~ {m/2)tc/U ~ 0.47^V*d- For instance, for the NV- 
diamond system (7^^ — 25 ns and — 32 ^jls 2M)^ 
we have e ^ 3 x 10^^. 

Secondly, detector dark counts on either round of the 
EO can lead to a spurious 'success' of the EO, which can 
reduce the fidelity of the entanglement. For existing APD 
detectors, dark count rates are typically Fdc < 500 s^^ 
[ifij i . Dark counts can be made negligible by observing the 
detector output only for the window i^ait ~ 3F~[^.^^ 1 
ns for NV-diamond) . The probability of a spurious count 
is therefore pdc — Tdc^wait ^ 10~^. Thus dark counts 
should have a negligible effect on the cluster fidelity. 

Finally, imperfect mode matching of the photons emit- 
ted by the matter qubit-cavity systems reduces the fi- 
delity, because the photons carry information regarding 
their origin. Non identical central frequencies, different 
polarizations, and spatio-temporal mode shapes of the 
photons can all reduce the fidelity. The frequency of 
the photons emitted from cavity i depends on the fre- 
quencies of both the ^ \e)i transition (w|e,i) and 
the cavity mode (a;cav,i)- The wxe/i's can be tuned in- 
dependently, e.g. by using local electric and magnetic 
fields to induce Stark and Zeeman shifts. The Wcav.i's 
can also be accurately and independently tuned, e.g. by 
using strain-tunable silica microcavities 1271 . or piezoelec- 
trically tuned fibre optic microcavities |28l |. The polar- 
ization of the emitted photons can be accurately matched 
using linear optical elements [i^. The spatio-temporal 
mode shapes of the emitted photons depend on the giS 
and KiS of the respective cavities. These parameters de- 
pend on the structure of the cavities, and hence are more 
difficult to calibrate once the cavities have been fabri- 
cated. However, we calculated that the EO is rather 
robust to mismatches in the giS and HiS (see Figl^D): 
mismatches of a few percent lead to a reduction in fidelity 
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FIG. 2: (a) Probability of success p versus tlie detector effi- 
ciency r] after both rounds of the entangling procedure. Each 
curve is plotted for a different value of the spontaneous emis- 
sion rate 7 = {0, 0.lFsiow, 0.2rsiow} (top to bottom). The 
other parameters used were g = 0.3 and ti = 1. (b) Fidelity 
of the entangled qubit pair for imperfectly mode-matched 
photons, calculated using the quantum trajectories formal- 
ism. Upper curve; fidelity for non identical leakage rates, 
X = ki/k2, taking gi = g2- Lower curveifidelity for non iden- 
tical coupling parameters, x — 51/(72, taking 52 ~ 0.3 and 

Ki = K2 — 1. 

of less than 10~'^. Deterministic sources of indistinguish- 
able photons, which have similar requirements to those 
needed for our scheme, are currently being developed by 
a number of groups [H US HI HI iHl . 

The next step towards scalable quantum computers is 
linking qubits together into cluster states, using the EO 
described above. A cluster state of qubits {qi, q2, ..qN} 
can be represented graphically by a collection of qubit 
nodes connected by edges connecting neighboring qubits, 
as depicted in Fig. A linear cluster of N qubits (a 
chain) may be represented in the form \C)i,,,n = (|T)i + 
U)i^2)(|T)2 + U)2^3).-(|T)w + U)w), where represents 
the Pauli phase- flip operation acting on qubit i. Such 
linear clusters can be grown using our EO, as we now 
describe. Given a cluster \C)2...n, qubit 1 can be added 
to the end of the cluster by first preparing qubit 1 in 
the state = -t- |t)i and then applying the EO to 
qubits 1 and 2. If the EO is successful, the resulting state 
isoftheform(|t)i|i)2±|i)i|T)2Z3)|C')3...Ar, depending on 
whether both clicks were observed in the same detector. 
This can be transformed into a cluster state by applying 
the local operations H1X2 or X1H1X2, conditional on the 
outcome of the EO (here Hi is the Hadamard operation, 
and Xi the Pauli operator implementing a bit flip). If 
the EO fails, the state of qubit 2 is, in general, unknown. 
However, measuring qubit 2 in the computational basis 
removes qubit 2 from the cluster, but projects qubits 
{3 . . . N} back into a pure cluster state. Therefore, failure 
of the EO causes the original cluster to shrink by 1 qubit. 

Repeatedly applying the procedure described above al- 
lows long chains to be grown. However, the theoretical 
upper limit on the success probability of our protocol is 
p = ^, and when the protocol fails, the chain shrinks 
by 1 qubit. Therefore, with only the procedure described 
above one can not create large clusters efficiently. If recy- 
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FIG. 3: Cluster states with arbitrary single-qubit measure- 
ments can perform a universal quantum computation, a) The 
qubits (circles) are entangled with their nearest horizontal 
neighbor via EOs (depicted by lines between the qubits) , and 
gates between computational qubits are incorporated by the 
vertical lines, b) Linear clusters can be joined together by ap- 
plying the EO between qubits Ai and Bi, and subsequently 
performing single qubit operations and measurements. 



cling of the clusters after a failure |2l| is not performed, 
the average number of EOs required to create a cluster of 



m qubits is A^eo = Z]i=i P 



P 



l-ml 



(l-p"-2)/(l-p). 



A way around this problem is a 'divide and conquer' 
approach, in which short chains are grown inefficiently, 
and then joined to a longer cluster using the EO together 
with local operations. The EO can be used to join two 
clusters, as shown in Fig. ^p. The initial state of two 
chains, A and B of length N and m respectively, may 
be written (|C)|^}|T)ai+^aJC){^}U)aJ(|T>Bi|C){b} + 
^B2li)Bi|C){s}), where \C){a} {\C){b}) represents a lin- 
ear cluster state of qubits {A2 . . . A^q} [{B2 ■ ■ ■ Bm}). 
These chains can be joined together by first perform- 
ing the local operation Xa^ , then applying the EO be- 
tween qubits Ai and and measuring Bi in the basis 
\±)bx = IDbi ± IT)bi- If the EO was successful, the re- 
maining qubits are left in the state ±|C){a} IT) \C){b}^ 
ZA2\C){A}\i)AiZB^\C){B}, where the first sign depends 
on the outcome of the Bi measurement, and the sec- 
ond sign depends on the outcome of the EO. Applying 
a local operation to qubit Ai yields a cluster state, of 
length 7V-|-m- 1, of qubits {An , ■ ■ ■ Ai,B2, ■ ■ ■ Bm}- If 
the EO fails, qubit Ai must be measured in the com- 
putational basis, and the original cluster state shrinks 
by 1 qubit. Thus the average length of the new cluster is 
L=p(iV + m- l) + (l-p)(7V-l). In order that the clus- 
ter grows on average, we require L > which implies 
that length of the short chains should satisfy m > 1/p. 

Chains of fixed length m can be grown independently 
using the EO, either by sequentially adding single qubits 
to the end of a cluster, or by joining sub-chains to- 
gether. Growing these m-chains adds a constant over- 
head cost to the cluster generation process. For ex- 
ample, growing a 4-chain (without recycling) requires 
on average p^^ + p^^ + p^^ applications of the EO, 
and each attempt to join such a chain adds on aver- 
age 4p — 1 qubits to the large cluster, leading to a to- 



tal cost of C4 — [p 



P 



p-^ + l)/(4p - 1) EOs 



per qubit added to the large cluster. A 5-chain can 
be grown by joining two 3-chains together. Joining 



such 5-chains to a longer cluster leads to a total cost 
ofCs = (2p-3-t-2p-2-|-p~i-hl)/(5p-l) EOs per qubit. 
To minimize these costs, the collection and detection effi- 
ciencies should be maximized. For example, for p « 0.24 
(or 77 = 70% with 7 = 0), we require to = 5, and we find 
C5 — 775. A modest improvement in detector efficiency 
dramatically reduces the overhead cost: for 77 = 85% 
and 7 = 0, we find C4 = 73.4. Note that there may be 
more efficient schemes for growing linear clusters using 
our EO (e.g. employing recycling of small clusters 0|) 
which yield lower overhead costs. 

In order to build linear chains into two-dimensional 
cluster states capable of simulating arbitrary logic net- 
works, cross links between linear chains must be con- 
structed [iOl . Such a link can be created by first using the 
EO to create an / shaped cluster (see Fig|3Ia)) offline, for 
some fixed cost. Provided the arms of this /-cluster are 
sufficiently long, the EO can be used to join the /-cluster 
to a pair of linear clusters with a high probability, and 
therefore create a cross link between the clusters. This 
leads to a constant overhead cost per cross link added to 
the cluster. Note that other methods for creating two- 
dimensional clusters, e.g. using microclusters [i^ or re- 
dundant encoding [2l|, have also been proposed. 

Our proposal has a number of very desirable features 
with respect to practical implementations. Firstly, our 
scheme requires only a simple level structure and single- 
qubit operations. Secondly, photon loss does not re- 
duce the fidelity of the entangled state of the qubits, but 
merely adds to the constant overhead cost. Thirdly, ow- 
ing to the simplicity of the optical networks used in this 
scheme, mode matching should be relatively straightfor- 
ward. Fourthly, the scheme is inherently distributed: in- 
dividual qubit-cavity systems can be placed in distant 
labs, and connected by optical fibers. This means that 
our scheme lends itself naturally to distributed appli- 
cations, such as quantum repeaters ^34] and quantum 
cryptography [s^- Finally, many of the techniques de- 
scribed here have been demonstrated experimentally, and 
the system requirements needed to create high-fidelity 
cluster states do not seem prohibitively restrictive. 

While preparing this manuscript, we became aware of 
an alternative scheme that may also be used for gener- 
ating cluster states of matter qubits |36| . We thank Tim 
Spiller and Bill Munro for valuable discussions and care- 
ful reading of the manuscript. The authors are supported 
by the E.U. Nanomagiq and Ramboq projects. 
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